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Abstract 

Recursion formulae of the iV-particle partition function, the occupation numbers 
and its fluctuations are given using the single-particle partition function. Exact re- 
sults are presented for fermions and bosons in a common one-dimensional harmonic 
oscillator potential, for the three-dimensional harmonic oscillator approximations 
are tested. Applications to excited nuclei and Bose-Einstein condensation are dis- 
cussed. 
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1 Introduction and summary 



In thermodynamics ideal quantum gases are usually treated in the grand 
canonical ensemble since all applications which were important in the past, 
like the electron gas, phonons or photons, deal with huge particle numbers, 
where the thermodynamic limit is applicable. Only the experimental attempts 
of the last years to investigate finite Fermi and Bose systems and to describe 
them in terms of thermodynamics called for theoretical effort. Interesting fi- 
nite Fermi systems are for instance nuclei, which behave like a liquid drop and 
therefore can undergo a first order phase transition [1]. On the low excitation 
site of the caloric curve the nuclear systems might be well described as an 
ideal Fermi gas in a common harmonic oscillator potential. 
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Small Bose systems became available through the development of traps. Here 
the focus is on the Bose-Einstein condensation which for instance could be 
found investigating alkali atoms in magnetic traps [2-4]. Again the system can 
be well described as an ideal quantum gas contained in an external harmonic 
oscillator potential. 

Theses finite quantum systems are characterized by a constant particle number 
and may be represented by the canonical or the micro-canonical ensemble 
depending on whether heat exchange with the environment is possible or not. 

In the following we present three formulae for the partition function of the 
canonical ensemble, one being already known [5]. All of them need the single- 
particle partition function as input. In a one-dimensional harmonic oscillator 
the result can be further simplified for fermions and bosons, leading to the 
interesting result that both have the same specific heat. 

In practice the exact formulae for the partition functions can be evaluated only 
for small particle numbers, say iV<,20. Therefore, approximations are derived 
in the third section, which may be employed for larger N. They are tested 
against the exact result for small N. In the last section we also discuss the 
examples of small excited nuclei and Bose-Einstein condensation. 



2 Recursion formulae and generating functions 



Throughout this article we assume that the Hamilton operator is a single 
particle operator, typically the kinetic energy and a common potential or a 
mean field 



N 

#=EM")- (1) 

n=l 

The partition function Zn for the canonical ensemble of N particles can be 
recursively built starting with the single-particle partition function [5] 



Z nW) = ^ £ (±l) n+1 Z 1 (n(3) Z N -M , Z (f3) = 1,^=^, (2) 

where the upper sign in the sum stands for bosons, the lower sign for fermions. 

A second method uses a generating function to obtain Z N . Consider the grand 
canonical partition function Q((3, z) 
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Q([3, z) = £ z N Z N ((3) , z = e^ , (3) 

Af=0 

where z is the fugacity and \i the chemical potential. Multiplying Q(/3, z) with 
a second function Y(f3, z) 



Y((3,z) = l + Y: (±l) n z n Z 1 (n(3) (4) 

n=l 

and grouping the result with respect to powers of z yields a differential equa- 
tion 



(1-Y((3,z))dz = dQCM 

Q(/3,z) 1 J 

which can be integrated obtaining 
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In {Q(P, ^)} = E (±!) n+1 - ^ == W z) (6) 



n=l 72 



or 



Q(P,z)=e W {W(P,z)} . (7) 

Expanding this function, which expresses Q{j3,z) in terms of Z\, into powers 
of z gives the Zn{[3) as coefficients. This derivation is equivalent to first differ- 
entiating Q{j3, z) eq. (3) with respect to z, inserting (2) and then regrouping 
the sum [6]. 

There is even a third method which represents Zn(/3) explicitly in terms of Z\ 



Zn(P) = ±* ^(±1)"-™ N[v] J] MnP). (8) 

is=0 n€K[v] 

We explain this method with the help of Z^(f3) and table 1. The method starts 
with the dual representation of v running from zero to 2 N ~ 1 — 1. This string [v\ 
is filled with zeros at the end to get as many digits as particles, in our example 
three. Then this string is divided into k[u] substrings each substring closing 
with exactly one zero, in our example for v — 1 the string 100 is grouped 
into 10 and 0. K\y\ is the A;[^]-tuple which consists of the lengths of those 
substrings. The last factor involved is N[u] which is the product of all integers 
from 1 to N taken to the power given by the respective digit in the string [u] . 
Adding up the results of our example leads to 
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Table 1 

Example of the explicit method eq. (8) for N = 3. 



ZM = ^ (Z 3 M ± 3Z 1 (P)Z 1 (2P) + 2Zx(3/?)) . (9) 

The proof of (8) is accomplished by regrouping the contributions v with re- 
spect to [u], where the first group contains all terms where [v\ ends with 00, 
the second where [u] ends with 010, the third where [u] ends with 0110 and 
so on. These groups can be identified with the summands of the recursion 
formula (2). The coefficients of the product Y\Zi{n(3) in (8) satisfy a number 
of combinatorial relations which will be investigated elsewhere. Here we only 
note, that the sum of all coefficients amounts to in the fermionic case and 
to 1 in the bosonic case. 

Following the very same idea that led to (2) one can derive recursion formulae 
for the occupation numbers of single-particle energy eigenstates | n ) with 
energy e n [7] 



N n ■= (( a + a n )) = ]T (±l) m+1 exp{-m (3s n } Z N _ m (P) . (10) 

Z N{P) m=1 

Any single-particle density Pn(P) (e.g. the spatial density) is given by 
1 N 

Pn(P) = E (±!) n+1 PiW) Z x (np) Z N _ n (P) , (11) 

Zn{P) n= i 

which follows from equation (10) using the density p\ in one-particle space. 
For the fermion case the fluctuations are well known to be 

5N n := (( (a+ n a n ) 2 )) - (( a + a n )) 2 = (( a + a n )) (l - (( a + a n )))(12) 
for bosons the fluctuations are given with the help of 
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2.1 Exact solution for the one- dimensional harmonic oscillator 



In the case of a one- dimensional harmonic oscillator potential the iV-particle 
partition functions can be given explicitly, they are (see also [8] for bosons) 



z„ (/j ) = exp(- W )n 1 _ exp( _„ flM . 

which differ for bosons and fermions only in Eq 



(14) 



E (N)=N^- for bosons 



E (N)=N'- 



for fermions 



(15) 
(16) 



The mean energy turns out as 



N t 
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(17) 



It is interesting to realize that ideal Fermi and Bose gases contained in a 
one-dimensional harmonic oscillator have the same specific heat 
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which is due to the equidistant spacing of energy levels and is, for instance, 
not valid in the case of a one-dimensional box. 



2.2 Exact solution for the one- dimensional box 



Since the formulae (2), (6) and (8) for the iV-particle partition function rely 
on the knowledge of the single-particle partition function it is always valuable 
to have Z\ explicitly. Considering a particle with mass m in a one-dimensional 
box of length L the energy eigenvalues are 
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and x :— 8 — 

2m V L J 



will denote the dimensionless inverse temperature. Then the partition function 
for the one- dimensional box can be calculated as 



n=l 



W = Ee-'4(«AO-l) = /^-| ■ (20) 



In three dimensions it is just the third power. t? 3 is the third elliptic -^-function 
and K the complete elliptic integral [9] 



& 3 (u, q) = 1 + 2 ]T q n2 cos(2rax) , (21) 

n=0 

K(m)=J* /2 d6 (l -m sin 2 (0))~^ . 
The variable m is related to the dimensionless inverse temperature x by 



^ 2 (0, e" 

m = 



M Q > e - x ). 



Mu,q) = 2q* E ? n(n+1) cos(2(n+ l)u) . (22) 



n=0 



In order to calculate the mean energy the first derivative of Z\ is needed which 
reads 



H 7 " 

' 2 AT§(m) [E{m) - {I - m)K{m)\ , (23) 



dx y/2 
with E 1 the complete elliptic integral 



E(m)= r /2 d9 (l -m sin 2 (#))^ . (24) 



3 Approximations 



5. 1 Fermions 



The approximation rests on the fact that a three dimensional harmonic os- 
cillator can be viewed as three independent one dimensional oscillators, one 
for each spatial coordinate. Then the idea is to apply the Pauli principle for 
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excited states only to one of the oscillators, i.e. directions in space. The ap- 
proximation of the partition function becomes [10] 



2 JV 



Z N (/3) = exp(-/3E (N)) 



1 — exp(— Ptvuj) j 



(25) 



N 



X 



n — 

n=l 1 



exp(— n fthu) 



with E (N) being the correct ground state energy. The approximate mean 
energy 



E N (T) = E (N) + 2N 



coth 



- 1 



(26) 



n=l 



coth n 



has besides the ground state energy three parts, two for the two directions 
where the Pauli principle is not applied and one where it is applied. Therefore 
the specific heat is a sum of two contributions for distinguishable particles in 
a common one-dimensional oscillator field and one for fermions in a common 
one-dimensional oscillator field. Although the partition function is not com- 
plete in this approximation the resulting mean energy turns out to deviate 
from the correct result by less then 1%, see example fig. 1. Only for iV = 2 
and iV = 3 the deviation is larger. The high temperature limit is correct. 




T (Wk R ) 




2 I 
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Fig. 1. L.h.s: mean energy of fermions in a three-dimensional oscillator, 
the solid line displays the exact result, the dashed line the approximation 
(26). R.h.s.: relative deviation between approximation and exact result. 



3.2 Bosons 



Unfortunately an approximation like (25) cannot be derived for bosons, here 
we try to derive an appropriate power series. For the ideal Bose gas in a 
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three-dimensional oscillator potential we have 



Z i = 77^-^3 > with x : = e "^ /2 • ( 27 ) 

(l — x z y 

Obviously, Zn will be a rational function of x. A closer inspection shows that 
it has the structure 



X 3N l n 
ZN ^ = TIN 7TZ 2^3 PN ^ > Where =HPm X 2 " 1 (28) 

is an even polynomial with non-negative coefficients where for N > 4 the 
degree is 2L N = (3N 2 - 7N + 10). The first (N + 1) coefficients p ,p 1 , ...,p N 
are independent of N, for example, 



P N ( X ) = 1 + 3x 4 + ■ ■ ■ for N > 1 . (29) 
Moreover, the following identities hold (the primes denote derivatives): 

P N (l) = (N\f (30) 

p'M = In(n-i)(n\) 2 

P&(1) = — (27N 2 - 23N - 8)N(N - l)(iV!) 2 

J. _ 

P^(l) = i(27iV 4 - 42iV 3 - 9N 2 + 167V - 20) N(N - l)(iV!) 2 . 
8 

These identities may be proved by deriving an analogous recursion relation 
for the P/v from eq. (2). 

In order to approximate Zn we consider the function 

f N (x)=\n[\n(P N (x))} , (31) 

which is analytic in a neighborhood of the real positive axis and has a loga- 
rithmic singularity at x = due to (29). Since /w(x) — 41n(x) behaves rather 
smoothly in the physically relevant interval x G (0, 1] it appears sensible to 
use a Taylor series approximation 



f N (x) =:41n(x) + g N (x) 
K a 

«41n(z) + \i x - • ( 32 ) 

i/=0 ZA 
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The values of the u-th derivatives ffc\l),v = 0,...,K, and hence of g u , 
can be expressed in a straight forward way as functions of Pj^\l) and, using 
identities of the form (30), finally as known functions of N, albeit of increasing 
complexity. We choose K = 3 and K = 4 and obtain approximations Z N , E N 
and c^, where the latter is compared to exact results in figure 2. 




T (fiu/k R ) T ttWk R ) 



Fig. 2. Specific heat of bosons in a three-dimensional oscillator, the solid 
line displays the exact result, the dashed line the approximation (32) 
with K = 3, the dotted with K = 4. 

One sees that the approximation already for small K reflects the exact be- 
haviour qualitatively correct. In the high and low temperature regime it works 
also quantitatively well, but the maximum is described only roughly both in 
position and heights. For larger particle numbers more and more terms of the 
Taylor expansion (32) would have to be taken into account. 



4 Applications 

4-1 Fermions 

Approximation (26) enables us to evaluate the low temperature behavior of 
the excitation energy of nuclei consisting of for instance A = 100 nucleons and 
compare it to the approximation of the free Fermi gas commonly employed in 
nuclear physics. Figure 3 shows on the l.h.s. as a solid line the result obtained 
for the common harmonic oscillator with a frequency of 

Huj « 41 MeV , (33) 

which is obtained if one demands a correct root mean square radius [11]. The 
dashed line displays the respective dependence for the free Fermi gas with 
an inverse level density parameter of 10 MeV. One clearly sees that for very 
small excitations the temperature rises much faster with excitation energy if 
calculated in the oscillator shell model. Then the curve bends earlier towards 
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Fig. 3. L.h.s.: relation between excitation energy and temperature for 
the harmonic oscillator (solid line) and the free Fermi gas (dashed line). 
R.h.s.: experimentally determined relation, figure taken from [1]. 

the high temperature limit c = 3ks, meaning that the temperature now grows 
less with excitation energy. In practice this observation is still beyond exper- 
imental verification as can be anticipated from the right hand side, where 
an attempt to determine the complete caloric curve of the nuclear liquid-gas 
phase-transition is displayed [1]. 

4-2 Bosons 

The observation of Bose-Einstein condensation of very cold gases contained 
in magnetic traps has stimulated much theoretical work in this field during 
the last years (see for instance [12-15] and references therein). Since the in- 
vestigated gases are very dilute and the trap can be modeled by a harmonic 
oscillator potential these systems can be viewed as ideal Bose gases in common 
harmonic oscillator fields. 

Before going into details let us first remind some known results [16]. Rigorously 
speaking a phase transition towards the Bose-Einstein condensate does not 
happen in an external potential since the chemical potential is always non- 
zero for finite temperatures. Only for a vanishing potential (infinite volume 
limit) it is zero for all temperatures below T c . It is also true that the phase 
transition does only occur for dimensions d > 2. 

Dealing with finite systems one has to relax the rigorous point of view and look 
for criteria that one would like to identify a "smooth Bose-Einstein conden- 
sation". A first natural criterion is the maximum in the specific heat [17,14], 
which is a kink in the free three-dimensional case because the mean energy 
jumps at T c . Looking at figures 4 and 5 (upper right corner) one sees that 
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Fig. 4. One-dimensional harmonic oscillator: mean energy, specific heat, 
ground state occupation number and its fluctuation for bosons (solid 
line) and distinguishable particles (dashed line). 
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T (Wk B ) T (fico/k B ) 

Fig. 5. Three-dimensional harmonic oscillator: mean energy, specific 
heat, ground state occupation number and its fluctuation for bosons 
(solid line) and distinguishable particles (dashed line). 



the specific heat exhibits a clear maximum in the three-dimensional oscillator 
whereas it does not in the one-dimensional oscillator. Using the specific heat 
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as a criterion for Bose-Einstein condensation the property that it does only 
occur for dimensions larger than two is maintained. 

A second possible criterion would be the ground state occupation number N 
and its fluctuation 5N . One could argue that the ground state should be 
macroscopically occupied below a certain temperature T c [12], but as figures 
4 and 5 (lower left) show, that happens also with distinguishable particles 
(dashed line). The same is true if one would take the existence of a maximum 
in the fluctuations of the ground state occupation number as a criterion. Also 
here the mere existence happens for distinguishable particles too, see figures 4 
and 5 (lower right). One also realizes that in this respect there is no difference 
between the one- or three-dimensional case. 
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Fig. 6. Fluctuations of the ground state occupation number in a 
one-dimensional (upper) and a three-dimensional harmonic oscillator 
(lower); lines from bottom to top starting with N = 5 in steps of 5. 
L.h.s. bosons, r.h.s. distinguishable particles. 

The possible key point in order to distinguish between bosons and distinguish- 
able particles is the dependence of the magnitude of the fluctuations on the 
particle number. For bosons in a one-dimensional harmonic oscillator the fluc- 
tuations <5iVo scale with N (see also [13]), whereas for distinguishable particles 
they scale with \fN. This can also be seen in fig. 6 (upper figures). 

In an experimental situation one of course knows whether one is dealing with 
bosons or not. The experimental interest then is to bring as many bosons as 
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possible into the ground state. Comparing the one-dimensional (fig. 6 upper 
figures) and three-dimensional case (fig. 6 lower figures) one realizes that for 
the same particle number this happens at higher temperatures in one dimen- 
sion. 
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